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[1] This paper is devoted to the parameterization of the Mueller matrix of oceanic water 
measured by Voss and Fry [1984]. Their matrix was obtained as an average of many 
observations in different parts of the world ocean, including the Gulf of Mexico and the 
Pacific and Atlantic Oceans. The results presented here can be used as input data for the 
theoretical polarized radiative transfer studies in oceanic waters. INDEX TERMS: 4271 
Oceanography: General: Physical and chemical properties of seawater; 4275 Oceanography: General: Remote 
sensing and electromagnetic processes (0689); 4552 Oceanography: Physical: Ocean optics; KEYWORDS: light 
polarization, light scattering, phase matrix 
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1. Introduction 


[2] Light propagation in oceans is studied in the frame- 
work of the radiative transfer theory [Mobley, 1994]. For 
this, one needs to know local optical characteristics of 
oceanic water such as absorption and differential volume 
scattering coefficients [Shifrin, 1988]. The main problem of 
the theoretical modeling lays in the fact that both concen- 
tration of suspended particles and their geometrical param- 
eters such as size, shape, and structure are highly variable 
both in time and space domains. 

[3] The situation becomes even more difficult if one is 
interested in the exploration of the light polarization in 
ocean water. This problem, for instance, is of importance 
due to a high correlation of the degree of polarization with 
the concentration of suspended particles [/vanoff; 1975], 
which potentially has important practical applications. The 
model of spherical particles, which is used even in up-to- 
date studies in ocean optics [Stramski et al., 2001], is totally 
inapplicable if one has a goal to describe the Mueller matrix 
[Bohren and Huffman, 1983; Kokhanovsky, 2001] of oce- 
anic water. For instance, the element P22 of this matrix is 
equal to P;, for spherical polydispersions. This is not 
confirmed by measurements both in open ocean and coastal 
areas [Kadyshevich et al., 1974]. 

[4] Clearly, the theory of light scattering by nonspherical 
particles should be used in this case. However, both the 
complexity of the nonspherical scattering theory and the high 
variability of ocean water microstructure make the solution of 
the problem highly complicated. With this in mind, we 
propose here the parameterization of the phase matrix of 
oceanic waters, based on experimental measurements. The 
experimental results used were obtained by Voss and Fry 
[1984], using an electro-optic light scattering polarimeter. 
Measurements were done on samples from the Atlantic and 
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Pacific Oceans and Gulf of Mexico. Only averaged exper- 
imental data for the matrix, obtained in different locations and 
time intervals, are taken for the parameterization. 

[s] The parameterization obtained can be used, for exam- 
ple, in the combination with analytical results of the 
polarized light transfer theory [see, e.g., Zege and Chaikov- 
skaya, 1996, 1998, 1999]. 


2. Mueller Matrix 


[6] The Mueller matrix relates Stokes vectors of incident 
and scattered light. It has the following general form 
[Bohren and Huffman, 1983]: 


Pu Py Po Pig 

r Pa Pu P3 Pa 

P= l 
P3; P32 P33 P34 (1) 
Pay Par Py Pu 


where elements P33, Pia, P3, Pa, P31, P32, Päis Pa, P34, 
P43 for oceanic water take very small values [Voss and Fry, 
1984] and can be neglected for most studies. This, of 
course, simplifies theoretical analysis. However, these small 
values in the right upper and left lower 2 x 2 submatrices 
(and also elements P34, P43) should not be considered as 
exact zeroes in general. They could be of importance in the 
solution of selected ocean optics problems, especially 
focused on the chiroptical spectroscopy of oceanic waters 
[Kokhanovsky, 2003]. Also, there are experimental indica- 
tions that these elements can have quite large values in 
single cases [Kadyshevich et al., 1974]. This could be due to 
orientation, anisotropy, or chirality of particles. We, 
however, will neglect all these effects and consider the 
matrix of the form 


Pi, Pin 0 0 

5 | Pa Po 0 0 

= 0 0 P3 0 (2) 
0 0 O Pu 
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This form has been confirmed by experiments of Voss and 
Fry [1984]. They also found that P}z ~ P21 and P33 ~ P44. 
Then the normalized Mueller matrix p = P/P,, takes the 
following form: 


~ | Pi2 pro 0 0 
P=\"0 0 pw 0 |) @) 
0 0 0 pss 


We are interested in the modeling of the dependencies of 
elements p12, P22, and p33 on the scattering angle 6. The 
dependence P,,(8) was studied by many authors [see, e.g., 
Petzold, 1972; Gordon, 1974; Shifrin, 1988; Stramski et al., 
2001; Mehrtens and Martin, 2002] and therefore is not 
considered here. Note that the degree of polarization p(0) = 
—p (8) [Kokhanovsky, 2003]. 

[7] The elements of the normalized matrix p of oceanic 
water are given by simple analytical equations, if only 
molecular light scattering is accounted for (no particle 
scattering) [Shifrin, 1988], 


p” (90°)sin? 0 


m 0 4 
Piy( ) 1 + p” (90°)cos? 0’ ( ) 
i _ p”(90°)(1 + cos? 0) 
AC) 1 + p” (90° )cos? 6” (5) 
+ 2p” (90°)cos 0 
P33(8) (6) 
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where p”(90°) = —pj3(90°) and the subscript “m” under- 
lines the fact that equations (4)—(6) describe molecular light 
scattering only. Note that p’”"(90°) = 0.835 at the wavelength 
equal to 436 nm and the temperature 20°C [Shifrin, 1988]. It 
is easy to show that the parameter p”(90°) gives the degree 
of polarization (at 8 = 90°) of initially unpolarized light, 
which is due to the scattering process by an idealized 
sample of oceanic water free from particles. 
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Figure 1. The dependence of the degree of polarization 
p(8) = —p;2(®) on the scattering angle © according to 
equation (7) (line) and measurements (dots). 
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Figure 2. Same as in Figure 1 but for p22(0), given by 
equation (8). 


[s] Generally speaking, the account for particle scattering 
modifies equations (4)—(6). Our task is to find analytical 
equations similar to equations (4)—(6), but for natural 
oceanic waters. For this we use experimental data obtained 
by Voss and Fry [1984]. We show experimental data by 
symbols in Figures 1—3. Lines correspond to following 
equations: 


p(90°)sin* 0 


p(0) = T+ p(90°)cos? 0” (7) 
pat AOSD 
pa(0) 2p(90°)cos 0 (9) 
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Figure 3. Same as in Figure | but for p33(0), given by 
equation (9). 
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[9] They differ from equations (4)—(6) only due to a 
different value of the degree of polarization at 90° (p(90°) = 
0.66), which was taken from Voss and Fry [1984, Table 2]. 
Also, the parameter z = 8 — 09 was introduced to describe 
the shift of the minimum in the dependence p22(09) from the 
angle 6 = 90°(compare equations (5) and (8)). We found that 
8) = 0.25 for the data given by Voss and Fry [1984]. 

[10] First of all, note that equation (7) works well for 
natural oceanic waters. Thus it is enough to know the degree 
of polarization at 90° to characterize the degree of polar- 
ization in the range of scattering angles 20—150°. For 
smaller and larger angles, where the degree of polarization 
is below 15%, however, the discrepancy increases. 

[11] Figures 2 and 3 show that equations (8) and (9) can 
be used also for modeling functions p22(0), p33(8) in the 
backward hemisphere. In the forward hemisphere, however, 
equations (8) and (9) contradict measurements. Therefore 
they should be modified. This necessity arises from the fact 
that equations (4)—(6) describe the dipole scattering mech- 
anism only. In the forward scattering range, however, the 
contribution of multipole scattering (which is due to large 
oceanic particles) is not negligible [Shifrin, 1988]. 

[12] This contribution can be estimated, taking into 
account that Pr. + P33 œ% P1; (see results of measurements 
in Figures 2 and 3) for small angles. We approximate 
elements P11, P22, P33 for multipole scattering at small 
angles by the exponential function €exp(—a8) with free 
parameters a and €, which is a standard approximation as 
far as the analytical radiative transfer in oceanic water is 
concerned [Zege et al., 1991]. Then it follows instead of 
equations (8) and (9), 


p(90°)(1 + cos?(6 — 09)) + £exp(—ad) 
1 + p(90°)cos? (0 — 0) + €exp(—ad) ° 


px (0) = (10) 


2p(90°)cos 0 + Eexp(—a) 
1 + p(90°)cos? 6 + Eexp(—ad) ` 


P33(8) = (11) 


The first term in the nominator of equation (10) is due to the 
dipole scattering and the second one is due to the multipole 
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Figure 4. Same as in Figure | but for p22(0), given by 
equation (10). 
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Figure 5. Same as in Figure | but for p33(8), given by 
equation (11). 


scattering. The same applies to the nominator of equation 
(10). A similar interpretation can be given to dominators. 
[13] Parameters a and € in equations (10) and (11) were 
obtained by fitting procedure, using Voss and Fry [1984] 
data. The result is: a = 4 rad~', £ = 25.6. The comparison of 
equations (10) and (11) with experimental data at these 
values of a, € and p(90°) = 0.66, 8 = 0.25rad is given in 
Figures 4 and 5. We see that the correspondence of 
experimental data and fitting curves is excellent (with a 
minor exception in the range of angles 40—60° for p22). 
[14] This makes it possible to state that the angular 
dependence of the normalized Mueller matrix of oceanic 
water can be described in a good approximation, if only 4 
parameters (e.g., p(90°), 8, a and £) are known. These 
parameters have a clear physical sense as it is shown above. 
Note that 9) = a! in the case presented here. Alternatively, 
one can use the following set of parameters: p(90°), Q%, 
P2290), P33(80). This is a major conclusion of this work. 
[1s] The accuracy of proposed approximate equations (7), 
(10), and (11) is studied in Figure 6. One can see that the 
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Figure 6. Relative errors of approximation (7), (10), and 
(11) for elements p12, P22, P33- 
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error of the parameterization is typically in the range [—5%, 
5%]. It increases for scattering angles, where normalized 
matrix elements approach zero. We can state that the relative 
errors of approximations (7), (10), and (11) are generally 
smaller or comparable to the errors of experiments involved. 
This confirms the fact that they capture correctly the 
behavior of functions pı2(0), p22(8), and p33(8) in natural 
oceanic waters. 


3. Conclusion 


[16] In conclusion, we proposed here equations (7), (10), 
and (11) for the parameterization of the elements of the 
normalized Mueller matrix (see equation (3)) of oceanic 
water. They can be used for the theoretical modeling (both 
numerical and analytical) of the polarized radiative transfer 
in oceanic waters. Another important point is that they 
suggest convenient analytical forms for the description of 
experimental data. The detailed relation of four free param- 
eters (p(90°), 9%, a, €) to the microstructure of oceanic 
waters should be clarified in future research. 
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